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Abstract 

The formation of interface from an initial sharp interface in polydisperse A/B blends is studied 
using the external potential dynamic method. The present model is a nonlocal coupling model 
as we take into account the correlation between segments in a single chain. The correlation is 
approximately expressed by Debye function and the diffusion dynamics are based on the Rouse 
chain model. The chain length distribution is described by the continuous Schulz distribution. 
Our numerical calculation indicates that the broadening of interface with respect to time obeys 
a power law at early times, and the power law indexes are the same for both monodisperse and 
polydisperse blend. The power law index is larger than that in the local coupling model. However 
there is not a unified scaling form of the broadening of the interface width if only the interfacial 
width at equilibrium is taken into account as the characteristic length of the system, because the 
correlation makes an extra characteristic length in the system, and the polydispersity is related to 
this length. 
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I. INTRODUCTION 



Polymers are often incompatible. When they are mixed together, an interface between 
them will occur. Because many properties of the blends are finally determined by the 
thickness of the interface and the concentration profile of the polymer across the interface, 
the equilibrium properties and dynamics of the interface attracted increased attention and 
gained widespread studies. The broadening of initially sharp interface between two types of 
polymers is of considerable interest. 

Self-consistent field theory is a powerful method in the study of inhomogenous equilibrium 
systems [1, 2]. It can give detailed distribution of the monomer concentration across interface 
and specify other thermodynamic and structure quantities of interest. In order to study the 
kinetic pathway towards its equilibrium state. A method called "dynamic mean-field theory" 
(DMFT) [3, 4] or "dynamic density functional theory" (DDFT) [5, 6, 7] was developed. It is 
a generalized time-dependent Ginzburg-Landau theory for conserved order parameters. The 
volume fraction (concentration) is often chosen as the order parameter in polymer blends. 
In this theory, the free energy is calculated using the mean-field approximation, i.e., in each 
time step, the concentration and the potential should satisfy the self-consistent equations. 
This free energy is more exact than the Flory-Hugins type and the Ginzberg-Landau- Wilson 
type free energy which arc both a local free energy plus the square gradient term of the 
concentration. In DMFT or DDFT, the Onsager kinetic coefficients are crucial quantities. 
Various chain dynamics and interactions will result different forms of kinetic coefficients. 
A simple one is the local coupling form which neglects the non-local interactions. In fact, 
sometimes the kinetic coefficients are assumed to be constants. If the dynamics of polymer 
chains are described by Rouse model, a different form of kinetic coefficient will be obtained, 
and the coupling between concentrations can be approximately described by Dcbyc function 
in a homogenous homopolymer melt [7]. In the reptation regime, many researchers argued 
that the forms of the kinetic coefficient were the same as that for Rouse dynamics in a 
homogeneous polymer melt [8, 9, 10, 11], i.e., proportional to the Debye function. However, 
Maurits and Fraaije [7] found that the kinetic coefficients of the two different dynamics are 
not the same though do not differ very much. Wang and Shi [12] studied the interdiffusion 
process between incompatible polymers using the Flory-Huggins type free energy plus the 
square gradient term of the concentration. The growing of the interfacial width W obeys 
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the general scaling form W{t) oc t°'. It is found that the exponent a is smaller that 0.25 
for all the cases in their studies. As the square gradient approximation is only valid near 
the critical point, it is appropriate to obtain a more exact free energy which should also be 
valid far away from the critical point. The DMFT is the desirable method. Using DMFT 
with the local coupling approximation, Yeung and Shi [4] studied the dynamics of polymer 
interfaces. They found that the exponent a is about 0.25 at early times, and is independent 
of the Flory-Huggins parameter x and the chain length. The interfacial width saturates to 
its equilibrium thickness at long times. Many experiments also focused on the interdiffusion 
process of a system with an initially sharp interface [13, 14, 15, 16, 17]. These experiments 
found that the broadening of the interface obeys an power law at early times, and it needs a 
very long time to saturate to its equilibrium thickness. It was also found that the exponent 
a was in a range from 0.25 to 0.4, and decreases as the temperature is lowered [13]. One 
argued that the quick increase of the initial broadening is expected to be controlled by the 
fast single-chain mobility, while the long time saturation to the equilibrium thickness is due 
to the large-scale hydrodynamic flow [17]. 

The dynamics of polymer interfaces are often described by the time evolution of the 
concentration of the chains across the interfaces. In DMFT, this can be done by solving 
the time-dependent Ginzburg-Landau equation which the concentration satisfies. As this 
equation is usually solved numerically in real space, the Onsager coefficient is approximated 
as a local coupling form for the simplicity and saving time in calculation [3, 4, 5]. Since the 
concentration is a conserved quantity and it is linear rated to the auxiliary potential field 
[18, 19], we expect that the auxiliary potential field is also conserved. The time-dependent 
Ginzburg-Landau equation with respect to concentration can be transformed to the form 
in which the auxiliary potential field is treated as the variable. Then we can obtain the 
external potential dynamic (EPD) model, as the equation is for the auxiliary potential 
fields. The EPD model is first proposed by Maurits and Fraaije. There are some advantages 
of the EPD method compared to DMFT method [20]. It incorporates a non-local couphng 
Onsager coefficient corresponding to the Rouse dynamics. In the spectral space, this can 
be realized through a local Onsager coefficient. Also, it is proved that numerically the EPD 
equation converges more easily and faster. In the present work we adopt the EPD method. 

Polymer chains in the real system are always polydisperse, and polydispersity may play 
important roles in determining the properties of materials. Polydispersity enriches the phase 
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behavior at equilibrium [21] and decreases the free energy barrier of nucleation in meta-stable 
state of polymer blends [22] . The effort of polydispersity on the profile of the interface at 
equilibrium are studied by Predrickson and Sides [23] in polydisperse polymer blends. In 
the present work, we studied the dynamics of interfaces in polydisperse blends. 

This paper is organized as follows. In Sec. II we present the derivation of the EPD model 
based on the Rouse chain dynamics. The methods for solving the EPD equations are also 
talked. In Sec. Ill the main results and discussion are presented. In Sec. IV we summarizes 
our conclusions. 

II. EXTERNAL POTENTIAL DYNAMIC MODEL 

In the present work wc consider an incompressible polymer blend of type A and type B 
linear, flexible homopolymers in which both species have polydisperse chain lengths. These 
two species of chains are modeled as Gaussian chains. To study the interfacial problems, it 
is convenient to work in the canonical ensemble. Por simplicity the distributions of chain 
length for both species are described by the same prescribed function P{N), where N is 
the degrees of polymerization, thus they have the same number- and weight-average chain 
lengths. Also both A and B monomers are assumed to have the same monomeric volume 
Po ^ and Kuhn lengths b. In this section we first derived the time evolution equation the 
concentration satisfies based on the Rouse dynamics in the polydisperse A/B blend. Similar 
derivation in the monodisperse case can be found [7] . Then this equation of motion has been 
transformed to the EPD form. The EPD equations are then solved numerically. 

A. Rouse dynamics 

Suppose that a polymer chain of type A with length N is subjected to an external force 
/A[RAr(s)], where RAr(s) denotes the position of monomer s. In the Rouse regime, as the 
correlations due to internal forces relax faster than the coarse-grained collective dynamics. 
Thus the chain can be considered as drifting with a constant velocity[7, 24] 



where Dq is the diffusion coefficient of a monomer. It is important to found from Eq. (1) 
that the drift velocities of chains with different lengths are different, while the velocities are 
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(1) 



the same for all beads in the same chain. The microscopic density of A species is defined as 

PA{r)^Yl E / c?s5[r-Riv(«)]. (2) 

N=l i=l ■''^ 

Here, riyuv is the number of A chains for chain length A^. The density can be alternatively 
expressed in the continuous limit by the following form, 

PA{r) = UA / dNP{N) / 6[r- R^(s)] , (3) 
Jo Jo 

where ua is the total number of A chains. According to the equation of continuity and using 
the chain rule, we obtain. 



^Pa{v, t) = -UA / dNP{N) / [r - R^(s)] • (4) 

ot Jo Jo ot 



We then replace dIlN{s)/dt by Eq. (1), and the external force can be expressed by minus 
the gradient of the chemical potential. By some calculation we obtain, 

d f foo p(N] r'^ 

-pA{r,t)=DonAVr-J dr' dN^^ ds ds'5[r-R^r(s)]5[r'-R^(s')]Vr'/XA(r'), 

(5) 

where pA is the chemical potential for A species. The concentration or the volume fraction 
of A species is defined by 0a(i") = (pa)/po, where (...) denotes the ensemble average. Then 
we can obtain 

d — r 

-M^,t) ^ Dr,(f)AV ■ J dr'J^ dNP(N)NgD(r-r',N)Vr'PA(r'). (6) 

Here — Dq/N„, is the number-average chain length of A species, 0^ is the average 
volume fraction, 0^ = uaNji/ PoV , where V is the volume of the system. The following 
relation has been used in deriving Eq. (6), ds ds'5[r — Rjv(s)]5[r' — RAr(s')] = N'^goir — 
r', N)/V, where ^'^(r — r', N) is the Debye function with chain length A^, and it is valid for 
Gaussian chains [24]. Similarly we can obtain the time evolution equation for 0s(r, t). 
However, (/)yi(r, t) and ^^(r, t) are not independent, they should satisfy the incompressible 
condition, 0^(r, + 0B(r, = 1. The incompressible condition can be incorporated to the 
dynamic equation through introducing a potential U (r) which can be added to the chemical 
potential. After eliminating this potential from the time evolution equations for 0a (r, t) and 
0B(r, t), we can obtain, 

-0(r, t)^V- j dr'A(r, v')y,,p^{v'). (7) 



Here, 0(r) = 0a(i") — is the concentration difference, A(r, r') = 

2Dn(t>A(t>B !o° '^^P{^)^9d{^ — is the Onsager kinetic coefficient. 0^ is the av- 

erage volume fraction of B species, 0^ = nsNn/ poV , where ub is the total number of B 
chains, /^^(r) = /U^(r) — Pb{y) is the chemical potential difference. In the self-consistent 
theory, the chemical potentials are given by 6F/pQd4>A and dF/po54>B- F is the free energy 
in the canonical ensemble: 

drx0A(r)0B(r) - po E / c?ra;«(r)0„(r) - ^ ^« / c^iVi'liV) lng„(a;„, iV). 

(8) 

Here, a;^ are the auxiliary fields conjugated to Qa[^a-i are the single chain partition 
functions of the chain length N for a species. 



B. External potential dynamic equation 

We have mentioned that concentration is a conserved quantity and it is linear rated to 
the auxiliary potential field, then the external potential field can be treated as a conserved 
variable. The dynamics of this potential can be described by the equation of the time- 
dependent Ginzburg-Landau form. The free energy of the system can also be written as a 
functional of the external potentials [2, 23, 25], which has the form of 

^ jdv -^^^^j dNP{N)\nQA[pi^-H-,N] 

roo 

- UB dNP{N)lnQB[l^+ + l^-,N]. (9) 
Jo 

Here, //+ = {cua + <^b)/2, = {cub — <^a)/2 are treated as the external potentials. The 
external field p,- couples the the concentration difference (pA — <Pb, while /x+ couples to the 
total concentration + As there are two external fields in H, the numerical evaluation 
of the functional integral on the fields is difficult, often a saddle point approximation will be 
made when the fields p+ is integrated [19], which means SH[p,^, p_]/Sp+\^^^^* — 0. This 
is consistent with the incompressible condition, i.e. 

cj)A[p*+ - P-] + + 1^-] = ^- (10) 

The external potential dynamics can be expressed as 

|,.W^V./dr'A,p,(r,r')V„^^fal. (11) 



Here, Aepd(i", r') is the Onsager coefficient in EPD model, and it is related to A(r, r'). The 
exphcit expression of A(r, r') can be obtained from Eq. (7). According to the chain rule 

dcf>{r,t) r , 5</.(r,t) a/x_(r',t) 



f d<p[r,t) Oii_[r,t) 

J'^Unir'A) dt ' ^^^^ 



dt J 6fjL-{r\t) 

where 5//+ /dt — has been used. According to the hnear response theory, using the saddle 
point approximation, for Gaussian chains at any time we can obtain 

Po Jo 0UA{r) 

and 

'^'^^(r) ^ _nA r^j^rprm SHuQa ^ riA 
SujA{r') Po Jo SujA{r)SujA{r') poV 

Similarly, we can obtain the expressions for 5(f)B{r)/5u;B{r'), and 6(f)A{r) / 6ujB{r') = 
5(f)B{r) / 5u;A{r') — 0. For the incompressible Gaussian chain system, after some calcula- 
tions, we can get 

*^(r) _ _, y^)-^B(r)] ^ ^^^^^ ,NPm^,M. - r'. N). (15) 



-'■^ rmPlN U = rdiVPWiVV(r-r'.iV). (14) 

Po -^0 dujA{r)oujA{r') PqV Jo 



5p4r') 5[u;A{r')-u;B{r')] ^^^^ io ' ' N, 

Insert Eq. (15) into Eq. (12) and compare with Eq. (7), and also through a Fourier trans- 
formation, the EPD equation can be obtained in the spectral space 



dp_{q) D^x /o°° dNP{N)Ngn{ci, N) 



dt 2 J,->dNPiN)^j^{ci,Nf 

where (pA and (pB are functionals of p\ and //_, and fi*^_ is determined by the incompressible 

condition, i.e., Eq. (10). The Onsager coefficient in spectral space has the form of 

, , , D,.x J,rdNP{N)Ngo(q,N) 
''-°<''> - —S-dNP(N)SLjMci,Ny 

where goifl, — 2(x-|-e~^ — l)/x'^, x — RgQ^, and Rg — Nb'^/6 is the radius of gyration. In 
the monodisperse hmit, P{N) = S{N — Nn), then AepdIq) = Dox/2Nn, which is consistent 
with the result obtained by Miiller and Schmid [19]. In the present work, we adopt the 
continuous Schulz distribution for reflecting realistic chain length distribution. However, it 
is not intrinsic to our topics. The Schulz distribution has the form of 

P{N) = ^l^^^^ ' exp[-(fc + 1)N/N^], (18) 

where is the weight-average chain length. A; is a parameter related to the polydispersity 
index, a smaller k corresponds to a more polydisperse distribution, and the infinity of k 
corresponds to the monodisperse case. 



X 



(16) 



7 



C. Numerical calculation 



In the present work, we consider a one-dimensional system, and 0^ = 0g = 1/2. The 
EPD equation (Eq. (16)) and Eq. (10) are closed and solved numerically. We adopt the 
numerical method developed by Ceniceros and Predrickson [2, 25], which is a semi-imphcit 
Seidel relaxation scheme. In the present work we studied the healing process of an initially 
sharp interface. The initially sharp interface was prepared by a tangent function, and for A 
species, it is 



(^Ai{x) = h -TT—r tanh 

2 2tanh?7 



jycos 



27rx 



X 



(19) 



where (pAp and (pAm are the concentrations at the boundaries in a phase-separated system 
at equilibrium, p denotes the A-rich boundary, while m denotes the A-poor boundary, 
is the length of the system, x ranges from —lx/'2- to /a;/2. The parameter t] determines the 
sharpness of the interface, and a larger 77 corresponds to a sharper interface. In the present 
work, we choose 77 = 100 which is large enough in our study. The initial concentration for 
B species is (f)Bi{x) — 1 — (f)Ai{x). However, in order to proceed the evolution of the EPD 
equation, we need the initial external potentials /i/_(a;) and fii+{x) which target the initial 
concentrations. This is realized by numerically solving the following equations 

(pAlfJ-i+ix) - ///-(x)] - (f)Ai{x) = 0, 

Ml^i+i^) + - (t>Bi{x) = 0. (20) 

We also use the semi-implicit Seidel relaxation scheme. 

In the calculation of the concentrations 0Q(r) {a =A,B), the infinite integrations are 
performed using a Gauss-Laguerre quadrature formula. The concentration (see Eq. (13)) 
can be expressed as 

(21) 

where M = {k + l)iV/iV^. The Gauss-Laguerre quadrature formula is 

n„ 

/ dMe-''f{M)^J2\f{Mi), (22) 

where the abscissas (Mj) and weights (Aj) can be find in a mathematical handbook. This 
formula converges very rapidly, and we find that 8 points are sufficiently accurate for all 
cases. In the present work, we adopt — 8. 
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The relations between the partition functions of the single chain (Qa) and the end- 
integrated propagators (qa) are Qa[N,u!a\ — I dxqa{x, N)/lx, where Qa satisfy the modified 
diffusion equations, which have the form of 

= J— u;^{x)qa{x,s), (23) 

with the initial conditions qa{x,0) — 1. The modified diffusion equations are solved nu- 
merically using the pseudo-spectral method [2, 23] with periodic boundary conditions (This 
results two A-B interfaces). This method is unconditionally stable in any number of space 
dimension and have higher accuracy [26] than the Crank-Nicholson semi-implicit schemes in 
our experience. For the convenience some quantities are scaled in the calculation 

X x/Rg^, s ^ s/N^, t -> tDjRl^, (24) 

where Rgyj — Nyjb'^/6. We also used x^w instead of x to characterize the incompatibility. 
The length of the system were chosen as Ix = 40, and the number of spatial grid was 
Nx = 512, which means a mesh size of Ax = 0.078125. In solving the modified diffusion 
equations the prescribed contour steps were As < 0.0025, which ensured a eight-figure 
accuracy of the concentration profile. The time step was chosen smaller at early time and 
larger approaching equilibrium at later times. We also conformed that our results were 
independent of the mesh size and the system size by varying Ax and l^- 



III. RESULTS AND DISCUSSION 



In order to study the effect of polydispersity on the formation of interface, we first in- 
vestigated the evolution of interface in the case of monodisperse A/B blend. The model we 
used in the monodisperse case was a simplified one of the present polydisperse model, which 
is also a nonlocal coupling one. 

Figure 1 shows an example of the time evolution of the density profile in the monodisperse 
case, in which the volume fraction of A species (pAi^) at different times for xA^=2.5 is 
presented. The initial interface is very sharp, and it evolves to the equilibrium one eventually 
at long times. In this process, as show in Fig. 1 the 0-enhanced and -depleted bumps 
appear, which were also observed in the local couphng model [4]. This can be understood 
by symmetry as explained by Yeung and Shi [4]. If we consider the interface at x = — 10, and 



9 



the locations at x — —20 and x — correspond to the phase boundaries at infinitely far away. 
The chemical potential is expressed ets /i — 5F/5(f) (0 means 0a), and it equals to zero at 
X — —20 and x = 0; also the chemical potential is an odd function with respect to x — —10, 
hence = at a; = — 10. There must be extremum values of the chemical potential between 
the interface and the phase boundaries, which is a maximum value on the A-rich side of the 
interface and a minimum value on the A-poor side. The diffusion of the order parameter 
satisfies jx — —M((f))dijL/d(f), where jx is the current and M(0) > is a dynamical coefficient. 
Near the interface, (f) is transported from the 0-rich side of the interface to the 0-poor side. 
However, the current changes its sign where the chemical potential is extremized. Therefore 
at some place far away from the interface, (p must be transported away from the interface 
on 0-rich side and towards the interface on 0-poor side. This is why an enhancement of 
above the equilibrium value on the A-rich side and a depletion on the A-poor side appear in 
the calculation. The calculation indicates that bumps are larger at smaller while they 
are smaller at larger xN, and in the course of time the bumps move to the phase boundaries. 

What we are interested in is the growth law of interfacial width with respect to time. 
The density across the interface is not homogeneous, thus the definition of the interfacial 
width is not unique. As a quantitative measure of the interfacial width, we followed Yeung 
and Shi, and also Steiner et al. [13] and chose the inverse of the maximum slope of 4>a at 
the interface 



Here, A0a is the difference between the bulk value of 0^ in the two equilibrium phases. 
It is obvious that this definition is sensitive to the local structure of the interface. Other 
definition, e.g. taking into account the entire structure of the interface, can alternatively be 
chosen, however, it dose not change the scaling relations. 

Figure 2 shows the interfacial width defined in Eq. (25) as a function of time for different 
XN. The interfacial widthes are divided by their equilibrium ones ^. From the figure it 
can be seen that only at the very beginning the broadening of the interface obeys a power 
law with respect to time. In order to find this relation we rescale the time variable, which 
means a shift of the profile along i-axis. As shown in Fig. (3), if W{t) is sclaed by ^ and 
the time also scaled by ^, all the data from Fig. 2 approximately collapse onto one single 




(25) 
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master curve. It means that W{t) satisfies the following equation 



(26) 



where /mono is a dimensionless function with /mono(oo) = 1 and (3 = 1.0. It can be found 
from the figure that at the early times, W{t) satisfies the following equation W{t) oc ct°', 
where c, a time independent constant, is related to ^, and a is the power law index. Prom the 
data fitting we found that a approximately equals to 0.38. A direct scaling analysis using an 
expansion of the free energy to the second-order in V0 by the Cahn-Hilliard dynamics gives 
oc or oc t^^*. An more exact form of the free energy, the self-consistent mean 

field free energy with the Cahn-Hilliard dynamics also demonstrated W oc t^^"^. It is easy to 
understand this relation since there is only one characteristic length ^ the equilibrium width 
of the interface in the approximated Cahn-Hillard dynamic system. In the present model, 
we take into account the nonlocal coupling between segments of a single chain, hence, there 
exists another characteristic length, the correlation length. It is because of this correlation 
length, a simple scaling analysis with one characteristic length can not predict an exponent 
of 0.38. An prefactor of c is also a reflection of the existence of another characteristic length. 
This will be more precise in the polydisperse case, which will be talked later. An exponent 
of 0.38 lager than 0.25 means that the correlation between segments quicken the broadening 
of the interface. 

In the following we consider the formation of interface in a polydisperse A/B blend. In 
the process of the evolution of the volume fraction of A species from its initially sharp 
distribution to the equilibrium diffusive distribution, the enhancement of 0^ above the 
equilibrium value on the A-rich side of the interface and a depletion on the A-poor side are 
also observed in the numerical calculation. The magnitude of these bumps becomes smaller 
when xN'w becomes larger, and also the bumps move to the phase boundary as t grows. It 
is similar to the case that in monodisperse blend, and also can be understood by symmetry. 

Figure 4 shows the evolution of interface with respect to time in the polydisperse polymer 
blends, the interfacial width is divided by ^. The polydisperse parameters are chosen as A; = 
1, 2, and 3, corresponding to the case that the polydispersity is reducing. The monodisperse 
case corresponds to an infinity of k. The incompatibility parameters are chosen as xA^'w = 
3.5, 4, 6, 8 from weak segregation region to strong segregation region. Prom the figure it is 
seen that only at very early time there is a power law between the interfacial width and time. 
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In order to find this relation we also rescale the time variable. However, in the simulation we 
found that different forms of scaling of time were needed for different polydispersities. This 
is as expected since the extra characteristic lengths, the correlation lengths are different 
in blends with different polydispersities. There is not a unified scaling expression which 
contains only one characteristic length. Figure 5 shows the detailed information of forms 
of time scaling under which the data approximately collapse onto one single master cure 
for different polydispersities, respectively. Prom Fig. 5 (a) for the case oi k — 1 it can be 
found that if W/^ is plotted as function of t/^^-^, all data for different incompatibilities 
approximately collapse onto one single fine; from Fig. 5 (b) and Fig. 5 (c), one can see that 
the scaling forms of time arc t/^^'^ and t/$,^'^ for /c = 2 and k = 3, respectively. It means 
that for a polydisperse polymer blend the broadening of interface satisfies 

Wit)=Mt/e), (27) 

where fk is a dimensionless function with /fc(oo) = 1 and k denotes the polydispersity. From 
the simulation it is obvious that the parameter (3 decreases as the polydispersity decreases, 
and it approaches the value of monodispcrsc case 1.0 when k goes to infinity. Equation (27) 
demonstrates that there is not a unified scaling form of the broadening of the interface if we 
only take into account the equilibrium width of interface as the characteristic length. An 
unified scaling form which can include different case of polydispersities may be constructed 
if two characteristic lengths, the equilibrium width of interface and the correlation length, 
are taken into account. However, it is out of the present work. At the early time for all 
the polydisperse case the broadening of interface with respect to time satisfies W{t) oc t", 
where the power law a ^ 0.38 which is the same as that in the monodisperse case. Thus we 
can conclude that the nonlocal coupling quickens the diffusion of polymer chains across the 
interface, while the polydispersity has no effect on the power law index. 

IV. CONCLUSION 

The formation of interface from an initially sharp one in polydisperse A/B blend is 
studied using external potential dynamic method. The chains of Both A and B species are 
polydisperse, and their chain length distributions are described by the continuous Schulz 
distribution. Based on the Rouse chain model, we derived the time evolution equation the 
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concentration satisfies. In the Rouse dynamcis, the drift velocities of chains with different 
lengths are different, while the velocities are assumed to be the same for all beads in the 
same chain. The collective Rouse dynamics model is a nonlocal coupling model. This 
coupling is reflected by the Onsager coefficient which is related to the Debye function. The 
time evolution of the concentration then transformed to the EPD form, because of which 
are numerically more efficient to solve. The EPD equations are solved numerically using 
a semi-implicit Seidel relaxation scheme. Our calculation indicates that the broadening of 
interfacial width with respect to time obeys a power law at early times, and the power law 
indexes are the same for both monodisperse and polydisperse blend. This means that the 
polydispersity does not affect the form of this power law. The power law index is larger 
than that in the local coupling model. It is obvious that the nonlocal coupling quicken 
the diffusion of polymer chains across the interface. However there is not a unified scaling 
form of the broadening of the interfacial width if only the equilibrium interfacial width is 
taken into account as the characteristic length of the system, because there exists another 
characteristic length, the correlation length. The correlation lengths are different in blends 
with different polydispersities. An unified scaling form which can include different case of 
polydispcrsitics may be constructed if two characteristic lengths, the equilibrium width of 
interface and the correlation length, are taken into account. This will be studied in the 
future. 
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FIG. 1: Time evolution of the density profile in the monodisperse case for = 2.5. The (pA- 
enhanced and -depleted bumps appear, which can be explained by symmetry. 
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FIG. 2: The broadening of interface with respect to time for different x-^ in the monodisperse case. 

The inter facial width is scaled by ^. It shows that at early times the interfacial widthes broaden 
fast, and they saturate to the equilibrium ones at very long times. 
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FIG. 3: Data from Fig. 2 plotted in a scaled form. Here the interfacial width is scaled by ^, the 
time is also scaled by ^. The data approximately collapse onto a single master curve. At early 
times, the broadening of interfacial width obeys a power law, and the power law index is about 
0.38. 
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FIG. 4: The evolution of interface with respect to time in the polydisperse case for different x^x 
and different polydispersities. The interfacial width is scaled by ^ 
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FIG. 5: Data from Fig. 3 plotted in a scaled form. Here the interfacial width is scaled by ^, the 
time is also scaled by for A; = 1 (a), by ^^'^ for = 2 (b), and by ^^'^ for A; = 3 (c). The 
exponent of ^ approaches the value that in the monodisperse case as k goes to infinity. The data 
approximately collapse onto a single master curve for different polydisperse case, respectively. At 
early times, the broadening of interfacial width obeys a power law, and the power law index is 
about 0.38 which is independent of the polydispersity. 
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